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Abstract 

To describe ultrarelativistic heavy-ion collisions we construct a three-fluid hydro- 
dynamical model. In contrast to one-fluid hydrodynamics, it accounts for the finite 
stopping power of nuclear matter, i.e. for nonequilibrium effects in the early stage of 
the reaction. Within this model, we study baryon dynamics in the BNL-AGS energy 
range. For the system Au+Au we find that kinetic equilibrium between projectile and 
target nucleons is established only after a time t e ^j M « 5 fm/c ~ 2Ra u /icm- Ob- 
servables which are sensitive to the early stage of the collision (like e.g. nucleon flow) 
therefore differ considerably from those calculated in the one-fluid model. 
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1 Introduction 



One of the central interests of heavy-ion physics is the study of the equation of state (EoS) of 
hot and highly compressed nuclear matter [0 . The success of the hydrodynamical model 0, [| 
in describing the gross properties of heavy ion reactions at BEVALAC energies indicated that 
a collective and equilibrated state of highly compressed nuclear matter is indeed created in 
collisions of heavy nuclei. In particular, experiments (qualitatively) confirmed H 

1. enhanced particle emission into preferred solid angles in collisions of a heavy target 
penetrated by a light projectile with velocity larger than the speed of sound in nuclear 
matter, indicating the occurrence of a mach cone || ||, 

2. collective flow of matter in the reaction plane [bounce-off) and perpendicular to the 
reaction plane {squeeze-out) in semiperipheral collisions of heavy nuclei || 

3. the linear increase of the pion multiplicity with beam energy due to entropy production. 
In this model entropy is produced in the compression shock front, the subsequent 
expansion proceeds via simple waves which - in the ideal fluid limit - conserve entropy 
and thus the pion multiplicity ||. 

These successes motivated the application of the (relativistic) hydrodynamical model to 
higher-energy heavy-ion reactions. At very high baryon densities and / or temperatures 
a phase transition from ordinary hadronic matter to a quark-gluon plasma (QGP) is ex- 
pected |], ||. In the (one-fluid) hydrodynamical model the energy and baryon densities 



necessary for this phase transition are already reached in the BNL-AGS energy regime [[LC 
Within this model, the effects of a (first order-like) phase transition lead to 

1. a plateau in the excitation function of (pt) [|TT1 , 



2. a plateau in the excitation function of the pion multiplicity [TT2 

3. a local minimum in the excitation function of the collective nucleon flow (pf r /N) |~3] 
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4. a prolonged lifetime of the system at the "softest point" of the EoS |14 . 



However, the question arises whether this model is still valid at ultrarelativistic energies 
(AGS and above), where these effects are predicted. In the one-fluid model instantaneous 
thermalization of projectile, target, and produced particles is assumed. Therefore, momen- 
tum conservation leads to immediate complete stopping at the contact surface (the nuclei 
do not penetrate each other) and then energy conservation requires that all the kinetic en- 
ergy is converted into internal excitation of the stopped matter. As a consequence, huge 
(midrapidity) energy densities, temperatures, and baryon densities are produced. Thus, the 
conditions for the phase transition are already reached at comparatively low bombarding 
energy. 

In this paper we extend the hydrodynamical model in such a way that thermalization 
of projectile and target does not happen instantaneously. Within our model, projectile and 
target decelerate gradually while penetrating each other. We study the influence of these 
nonequilibrium effects on baryon stopping, directed baryon flow, and the compression of the 
nuclei. Based on these results, we comment on the validity of the one-fluid hydrodynamical 
model at ultrarelativistic energies and our expectations concerning the above mentioned 
predictions of this model. 
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2 The three-fluid model 

In this section we present the physical motivation for introducing three fluids and discuss in 



detail the interaction sourceterms (see also refs. [15, 16, PJJ). 



2.1 Physical motivation 

The description of the early stage of a heavy-ion collision within the framework of ideal 
one-fluid hydrodynamics tends to become more and more unrealistic, when nuclei collide at 
AGS-energies (y/s ~ 5AGeV) and above. The major reason for that is that the assumption 
of instantaneous local thermalization of projectile, target and produced particles fails at 
ultrarelativistic bombarding energies. This is due to the finite stopping power of nuclear 
matter. As a consequence, one-fluid hydrodynamics overestimates the baryon and energy 
densities at midrapidity during the early stage of the reaction. This in turn yields a high 
pressure in the central region, which produces a large sideward flow in the reaction plane 
(the so-called bounce-off) |I0| . 

To get a rough estimate of the compression at midrapidity, one may employ the one- 
dimensional shock model ]TB| with an ultrarelativistic ideal gas EoS (p = e/3) for the shocked 
matter. The baryon and energy densities are then given by 

n/n = 47 cm - 3/-f cm , (1) 

n ~ 0.16 fm~ 3 and e — 0.15 GeV/fm 3 denote ground state baryon and energy density of 
nuclear matter. The linear increase of the central baryon density with 7 cm oc ^ becomes 
questionable for 7 cm 3> 1, since the assumption of instantaneous thermalization of projectile 
and target matter ceases to be valid. 
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To account for nonequilibrium effects in the early stage of the reaction we introduce 
several fluids in our model. Considering collisions of single protons at that high energies, 
one observes a strongly forward-backward-peaked rapidity distribution Each proton 

loses about one unit of rapidity and is deflected only by a small angle. The projectile and 
the target protons remain separated in phase space, if the incident energy is high enough, 
and can be identified as the projectile respectively target proton even after the collision. 
This means that to first approximation the two baryon flows (projectile and target) do not 
exchange baryon charge. For this reason we consider two fluids that correspond to projectile 
(fluid 1) and target nucleons (fluid 2). Due to inelastic nucleon-nucleon scattering, energy in 
form of newly produced particles is deposited at midrapidity. Since the two nucleonic fluids 
do not populate midrapidity in the beginning of the reaction, if the initial rapidity gap was 
large enough, the newly created particles are also separated in phase space. This motivates 
to collect them in a third fluid - the fireball. 

In our approach projectile and target fluids penetrate each other and decelerate gradually, 
not instantaneously. Due to inelastic NN scattering the fireball is dynamically produced 
during the compressional phase of the reaction. Since the time scale of energy deposition 
is rather short, we neglect the rescattering of the nucleonic fluids with the fireball during 
this early stage. Later, when all three fluids populate the same phase space region and 
thermalize, they should be merged into one. Then, hydro dynamical expansion becomes the 
dominant process. 
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2.2 Multi-fluid relativistic ideal hydrodynamics 

The basic equations of relativistic one-fluid hydrodynamics are the conservation of energy 
and momentum, 

dfjT^ = . (2) 

If we assume an ideal fluid (i.e., neglecting dissipative effects), the energy- momentum tensor 
T^ v reads § 

= (e + p)u^u v - pg^ , (3) 

where e(x) is the energy density in the local rest frame, p(x) the pressure, u^(x) = 7(1, v) the 
4-velocity field (normalized to unity, u^u^ = 1), and g^ v = diag(l, — 1, — 1, — 1) the metric 
tensor. 

In addition, there may exist several conserved currents j'f , for which the continuity equa- 
tions 

W = (4) 
hold. We will consider only one such conserved current, namely the net baryon current 

f = nu^ , (5) 

where n(x) denotes the net baryon density in the local rest frame. 

As computational frame we choose the frame where the longitudinal velocities of projectile 
and target nuclei are of equal magnitude. In this frame equations (|2|||) read 

d t E + d % {Ev l ) = -difa*) , (6) 
d t M j + d i {M j v i ) = -djp , 
d t R + di(Rv l ) = . 
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E, M, R denote energy, momentum, and baryon density in the computational frame. They 
are related to the local rest frame quantities e, p and n via 

E = T 00 = 1 2 (e + p)-p , (7) 
M i = T i0 = 1 2 ( e + p)v i , 
R = j°=jn . 

The system (||J^) of coupled differential equations is closed via the assumption of local 
thermodynamic equilibrium and by specifying an EoS in the form p(e, n). 

In hydrodynamics for N fluids, we split the energy-momentum tensor and the baryon 
current of the total system into a sum of N terms, 

= rr + T^ + .-.+T/J" , (8) 

The total system is now subdivided into iV individual fluids. Note that eqs. (H) hold lo- 
cally, i.e. different fluids may coexist at the same space-time point. Each individual energy- 
momentum tensor T^ v and baryon current jf does not need to be conserved, since the iV 
subsystems may exchange energy, momentum and baryon charge: 

W = *T > ( 9 ) 

d,j? = S t (l = l,;N) . 

Fi denote energy and momentum sources of fluid I, and Si the baryon charge exchange for 
each fluid. Due to eqs. (@,|]), i.e. the conservation of total energy, momentum and baryon 
charge, we have the additional equations 

F? + FZ + .-- + F£ = , (10) 
S 1 + S i + -" + S N = . 



Following the discussion in section |2.1| , we choose N = 3 and identify fluid 1 with the 
projectile and fluid 2 with the target nucleons. Fluid 3 collects the energy loss of projectile 
and target nucleons due to inelastic collisions, and will thus be identified with the produced 
particles. At time t — (in the computational frame), the third fluid is not yet existent: 
j£(t = 0,£) = 0, Tf(t = 0,x) = 0. 



2.3 Coupling terms between the fluids 

We now discuss the couplings between the three fluids, based on the assumption that they 
are well separated in rapidity. 

We first discuss the interaction between projectile and target. The total one-particle 
distribution function can be written as the sum of two individual distribution functions for 
projectile and target nucleons, respectively: 

f(t,x,p) = f P (t,x,p) + f T (t,x,p) . (11) 

Due to the strongly forward-backward peaked cross section (see discussion in section 
2~T|) in high energy pp-collisions it is possible to identify the projectile respectively target 



nucleon even after the reaction. The nucleon scattered into the forward hemisphere belongs 
by definition to the same flow as the one previously heading in this direction. Consequently, 
projectile and target fluid do not exchange baryon charge: Si = S 2 = 0. From eq. ([[(]) and 
the fact that the fireball is also separated in phase space from the nucleonic fluids, it follows 
that the fireball remains net baryon free, S3 = 0. 

The source terms of the energy-momentum tensor Ff in eq. (||) describe the energy 
respectively momentum loss of fluid I per volume and per unit time. We define 

dV 



Fi = Tlx n 2 I V M0 Uer J d(TNN->NX ij?' -V)j : d(T N N^NX = ^NN^NX ' 



E' 
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pj,>0 



This is simply the collision rate per volume n\ n 2 VM0iier da for scattering into the invariant 
phase space volume d 3 p'/E' times the average 4- momentum loss (jp' — p) v in a single binary 
iViV-collision. <r^^ NX denotes the invariant differential cross section. p' u , p v are the 4- 
momenta of the particle after and before the collision, respectively. The integration extends 
over the forward hemisphere of the cross section only (pj, > 0). Quantities with a tilde refer 
to the center-of-mass frame of the binary NN collision. The n« denote the local rest frame 
baryon densities of projectile and target, as defined in eqs. (|) and (^). vutiier denotes the 



invariant relative velocity of the two fluids pO 



VM0ller = V («1 U 2li ) -1 . (13) 

The (-)-averaging over the individual particle distribution functions provides a smearing of 
the sourceterms (accounting for the thermal and Fermi momenta of the particles), especially 
of the M0ller velocity. It is negligible during the compressional stage of the reaction since 
the internal velocities within each fluid element are small compared to the M0ller velocity of 
the colliding fluids. We therefore ignore this averaging. Instead, when thermal and relative 
velocities of the colliding fluids become comparable in the later stage of the collision, the 
fluids are unified (see section ^75|) . 



In principle it is always possible to split the sourceterms in a symmetric and an antisym- 
metric part with respect to the fluid indices (1 «-> 2): 

U H 1 J exchange J loss ) 

f) t^ v = — f v — f v (^A\ 

2 J exchange J loss i 

d u Tf = 2f? oss . 

The antisymmetric term Exchange describes the exchange of energy and momentum between 
projectile and target fluid, while fi oss denotes the loss of energy and momentum transferred 
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to the fireball. 

Following [15, 16, |17| we split the integral in eq. (|12|) into two moments, 



Fx = ^n 1 n 2 v M 0iier [{P2 -pi) u (Jp(s) - {p 2 + P\) u cr E (s)} , (15) 
F% = ^n 2 n 1 v M 0Uer[(Pi ~ p 2 ) u ap(s) - (pi+ P2) v cr E {s)} ■ 

Since we neglect the (-)-averaging, the particle momenta are simply given by p\ 2 = m^u\ 2 . 
The moments of the cross section <je and <jp are defined as 



CTe(s) = J dfX N N->NX ( 1 - ^ j , CTp(s) = J d(TNN^NX U - ^ j (16) 
p'>0 ^ ' p'>0 ^ 

and parametrize the mean energy respectively longitudinal momentum loss in a single nucleon- 
nucleon collision. E, p are energy and momentum of the nucleon before the collision, while 
E', p' denote the same quantity after the collision. Again, the integration is carried out only 
over the forward hemisphere of the cross section (p'u > 0). 

Comparing eqs. (|T4] ) and eqs. ( pT5|) we see that oe corresponds to the symmetric term and 



leads to the production of the fireball. Analogously, op is the antisymmetric term describing 
the exchange of energy and momentum between projectile and target fluids. Also, from 
this form of the coupling terms it immediately follows that the third fluid is produced at 



midrapidity: in the CMS p~[ = —p 2 and thus F 3 = 0. The hydrodynamical eqs. ( |10D and the 



specific form of the sourceterms, eq. (|15"D, can also be derived from the Boltzmann equation 



We employ the parametrizations of ap t E from the two-fluid model of refs. [15, 16, 17]. 
There the cross sections of free binary iVW collisions (NN — > N*X) have been used to 
evaluate the moments op^E- The incident baryons are assumed to be nucleons, whereas in 
the final state we sum over all baryons N*. The rescattering with the fireball (Nit — > X) is 
not considered, since the contribution to the total cross section outside the A-resonance is 
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small at high bombarding energies. Since we focus on the dynamics of the nucleons in this 
paper, we do not expect this to play an important role. 



2.4 Equations of state 

To close the coupled system of differential equations (H) we assume local thermodynamic 
equilibrium within each fluid (but not in the total system !) and an EoS for each fluid. 

For the target and projectile fluids we employ an EoS of an ideal nucleon gas plus com- 
pression energy: 

p = ({e-E c n)+p c . (17) 

E c is the compression energy per nucleon and p c the compressional pressure. Since the 
nucleon mass is large compared to the typical temperatures of the nucleons, we choose 
the nonrelativistic limit ( = 2/3. For the compression energy, we employ the ansatz 

E c = -p—(n-n ) 2 + m N + W , (18) 

i.e., E c increases linearly with the baryon density n for n ^> n Q « 0.16 fm~ 3 . We assume a 
ground-state (in-)compressibility k c = 300 MeV and a binding energy of Wq = —16 MeV. 
The compressional pressure p c and the velocity of sound cs at temperature T = are given 
by 

Pc = -^ = n^ = ^-(n*-nl) , (19) 
an 1 an Iouq 

2 dp n/n 



b de n/n -l + 9(m N + W )/k c 
Ansatz (p~8|) thus ensures cs(n,T = 0) < 1 as long as k c < 9(mjv + Wo) ~ 8.3 GeV. 

For the third fluid we employ an EoS of an ideal gas of n, n, p, uj mesons with a phase 



transition to QGP at T c = 160 MeV ||21|| . However, here we do not discuss the dynamics of 
the third fluid. 
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2.5 One-fluid transition 

In the later stage of the collision the nucleonic fluids stop. Their relative velocity (M0ller 
velocity) is then comparable to the internal thermal velocities. The two fluids are no longer 
separated in phase space, so that the main assumption for a two-fluid region does not hold 



anymore. Moreover, the coupling sourceterms, eqs. (15), cease to be valid, since they do 
not account for thermal smearing and vanish linearly with the M0ller velocity. Thus, two 
infinitesimal fluid elements with equal fluid velocities, Ui = u 2 , would not exchange energy 
and momentum, even if their temperatures were different. Because the particle velocities 
are thermally distributed around the fluid velocity, there should be a non-zero coupling 
that further drives the two fluids into common local thermal equilibrium. Since we do not 



account for this, the two fluids are merged, if the M0ller velocity (|13|) is comparable to 
the root-mean-square velocity in a nonrelativistic degenerate Fermi gas or a nonrelativistic 
Boltzmann gas: 

( Av ) 2 Boltzmann = max , (20) 

( Av ) Fermi = mft{2e F (ni) /m N } , (21) 



where 



'3 \ 2/3 

e F (m) = ( -ttSJ /2m N (22) 



denotes the Fermi energy. In contrast to 0, y, [15f] , the transition is done abruptly without 
a smooth transition from two to one fluid. 

Under the constraining assumption that both fluids flow with a common velocity v after 
unification, the set of eqs. (0) is solved for the local rest frame densities e and n by inserting 
the total densities E = E 1 + E 2 , R = Ri + R 2 , M = Mi + M 2 . By virtue of the EoS, the 
total pressure p = p(e,n) is also known. Additive quantities, i.e. the pressure p = pi + p 2 , 
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the baryon density n = n\ + n 2 and the energy density e = t\ + e 2 are then redistributed 
to the individual fluids according to the ratio of the baryon densities in the computational 
frame f2|, |: 



Pi = p p Z g ' Wi = n g Z g ' €i = 6 b + g ' *e{l,2} . (23) 
Since this procedure changes the individual fluid velocities V{ to the common velocity v, the 
propagated quantities Ri, E iy Mi have to be redetermined from rii, e« and pi according to 
equations (0) with the common velocity v and the corresponding 7 factor. Let Ri, Ei, Mi 
be the resulting quantities. Our special choice of the redistribution, eq. ( p3|) , implies that 
the individual baryon densities in the computational frame are not changed: Ri = Ri 

It can be easily proven that the total quantities Ex + E 2 = E — Ex + E 2 , R± + R% = 
II /»'; • /.' 2 . .\/ ; • M-> M M\ ■ M-> fulfill the one-fluid hydrodynamical eqs. (|). That 
makes it possible to propagate the fluids separately. Even if that has no physical meaning 
for the individual fluid, one can tell what amounts in phase space originated from projectile 
or target, respectively. 

Moreover we will use this feature to simulate a one-fluid model with our three-fluid 
model, by switching off the production of the fireball (setting oe = 0) and forcing the local 
unification of the two nucleonic fluids at any instant. This is what we will call the one- 
fluid limit of the three-fluid model. Formally, the one-fluid limit corresponds to Av = 00 in 



eqs. (50,21 



In the present stage of model development, only the one-fluid transition of the nucleonic 
fluids is implemented. The unification of all three fluids is recommended and in preparation 
because much of the nucleonic fluids finally end up at midrapidity, where the third fluid is 
produced. This would slightly increase the pressure in the central region and enhance the 
directed and radial flow of the nucleons. 

14 



2.6 Technical realization of the model 

The hydrodynamical eqs. (§) are solved on an Eulerian grid using the SHASTA (SHarp 



And Smooth Transport Algorithm) fl23| , |24], ^5]. This algorithm was extensively tested in 



27 1 . Each of the three fluids is propagated on its own grid with eight different densities 



(R,M,E,n,e,p). For each fluid the three-dimensional propagation is performed using the 



method of time-step splitting (operator splitting) p6| , i.e. the 3-divergence operators on 
both sides of eqs. (@) are not summed over the contracting index i. One solves sequentially 
equations of the form: 

d t U + di(Uvi + f) = , toTi = x,y,z . (24) 

Here U stands for either R, Mi, or E and —dif denotes the corresponding right-hand side of 
eqs. (D), e.g. in case of U = E and i = x we have — d x f = —d x {pv x ). The z(beam)-direction 
is propagated first, those z-propagated densities are then propagated in ^-direction (impact 
parameter) and finally propagated in y-direction (out-of-plane) . This order is kept fixed for 
reasons of simplicity. By permutation the symmetry of the propagation could be enhanced 
26| . Obviously the propagational sourceterms —dif are included in the propagation. 

In contrast, we exclude in our calculations the coupling sourceterms (right-hand side of 
eqs. (U)) from the propagation. In our implementation the fluids are first propagated without 
coupling terms and their local rest frame densities are calculated according to eqs. (0). From 
these densities the coupling sourceterms (eqs. fljjP ) are calculated and added. Since this 
coupling changes the energy and momentum density of the individual fluids, eqs. (0) have 
to be solved again in order to accomplish the full time step. Applying this procedure it is 
easy to control the effect of the coupling, e.g. a cell cannot become more than empty by 
extracting energy from it. In such a case the coupling would not be applied to any of the 
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fluids. 

In order to minimize the Lorentz contraction we choose the equal-velocity frame as the 
computational frame. In the equal-velocity frame both projectile and target nucleus have 
the same velocity (with opposite sign). In the case of two nuclei of the same mass this is 
identical to the center-of-mass frame. Because of the occurrence of acausalities one should 



not use a different grid spacing in different spatial directions ||28|| . Therefore, we use an 
equidistant grid with grid spacing Ax = Ay = Az = 0.15 fm. This choice was found 
sufficient to provide convergence in the rapidity distributions. It ensures the correct number 
of NN collisions with the appropriate colliding energy y/s. Our time step is At = 0.4 Ax to 
fulfill the Courant-Friedrichs-Lewy criterion for the SHASTA. 



As a consequence of eqs. (0) the cell velocity is: 

M 



E + p 



(25) 



We call a cell acausal if E < 



M 



. This happens mainly because in the final step of the 
SHASTA the propagated density is distributed uniformly over a cell so that, on a too coarse 
grid, some part of the fluid is transported due to this smoothing procedure and not by real 
motion. Another possible origin is that due to the coupling of the different fluids too much 
energy is extracted from a fluid compared to the extracted momentum. We are correcting 
those acausal cells by rescaling the energy density E with the momentum density: 

E CO rrected = (l + KT 8 ) • |m| . (26) 

This corresponds to setting the corrected cells to a maximum velocity v max . By removing the 
acausal cells instead of correcting them as described above, we found for the Ax = 0.15/m 
grid used here that less than 5.0% of the initial baryon charge and 2.0% of the initial total 
energy was lost. 
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3 Previous models 



Before presenting results we briefly review some earlier work on which our model is based. 
The main difference to these approaches is that in none of the predecessor models the pro- 
duced particles were collected in a third fluid and propagated hydrodynamically. They were 
two- fluid models only (i.e., N = 2 in the above equations). 

To our knowledge, the first application of multi-fluid hydrodynamics to heavy-ion colli- 



sions was discussed in ref. ||22|| . Similar to our model, fluids 1 and 2 corresponded to the 
projectile and target nucleons, respectively. They could exchange energy and momentum 
but no baryons. The couplings between the two fluids were already derived using the under- 
lying concept of binary nucleon-nucleon interactions. However, only elastic nucleon-nucleon 
scattering (i.e., no particle production) was taken into account. In this model Ff + F$ = 
and thus the sum of T^ v and T% v was conserved. 

In ref. || a similar model with a more refined coupling (more appropriate for the upper 
BEVALAC energy range, -EfX — 2 AGeV) was presented. However, also in this model no 
midrapidity fluid was considered, i.e. Ff + F 2 M = 0. 



In contrast, the authors of [E0| investigated the case Ff + F^ ^ 0. Thus, the possibility 



of inelastic reactions in collisions between projectile and target nucleons was accounted for. 
However, the antisymmetric (with respect to exchange of the fluid indices 1 <-> 2) term in 
their couplings (diffusion term), 

F f -F£ = -2£>nm 2 «-ti£) , (27) 



involved a diffusion constant D, whereas we employ a -yi-dependent op, cf. eq. (^5j). More- 
over, in their model the symmetric term (drain term) was 

Ff + F 2 " = -2< 7 (n 1 < + n 2 <) , (28) 
17 



again with a ^/^-independent q. As already mentioned above, the radiative energy losses 
of fluids 1 and 2 were not collected in a third fluid and propagated hydrodynamically, in 
contrast to our model. Furthermore, the fluids were propagated in (1 + l)-dimensions only, 
thus calculations of the transverse flow of the nucleons were not possible. Finally, since 
merging of the two fluids into one fluid (when they reach thermal equilibrium) was also 
not considered, the expansion of the stopped and excited projectile and target fluids and 
collective flow of spectator matter could not be investigated. 

In the two-fluid model of [ 15| , [16 , |17| , the sourceterms were calculated using experimental 
data for the elastic and inelastic NN cross sections and do not contain free parameters like 
D and q in eqs. ( |2"T| ) and (|2"8|). There was no energy loss to the fireball below the pion 
production threshold and the energy deposition in the fireball increases with the cm. energy 
t/s of the colliding nucleonic fluids. This ^/s dependence of the sourceterms accounts for 
the multiple production of particles, mainly pions. However, the authors of [I5L [HI W7\ 
considered only two fluids: the energy loss due to the creation of the fireball was accounted 
for in the propagation of projectile/target fluids, but the fireball itself was not propagated 
hydrodynamically. In |JTJ the opposite point of view was taken by studying the production 
and hydrodynamic expansion of the third fluid in 5* + 5* at 200 AGeV (assuming some 
global interaction between the nucleons of fluids one and two), however, disregarding the 
compression and heating (and thus the hydrodynamic evolution) of the projectile and target 
fluids. 

In contrast, our three-fluid model propagates all three fluids hydrodynamically and treats 
the interactions locally in space-time. It thus allows to study spectra of produced particles 
like thermal pions or photons []2T| as well as the dynamics of the nucleons. In this paper we 
focus on the latter and thus neglect the scattering of the fluid of produced particles with the 
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projectile and target fluids. A more refined investigation of various particle spectra, as in 
32fl , is in preparation. 



4 Results 

In this section we present results on compression of projectile respectively target, baryon 
stopping and flow. We emphasize the deviations from the one-fluid limit which is achieved 
by setting <Je = (i.e. no energy loss due to particle production) and by merging the 
projectile and target fluids at each space-time point from the very beginning. 

A remark about the freeze-out should be added here. Usually a hydrodynamical treat- 
ment becomes unrealistic when reaching some local freeze-out criterion, i.e. for example the 
local temperature drops below the freeze-out temperature |£J or the local density becomes 
smaller than a certain freeze-out density. Fluid elements would then freeze-out locally (e.g. in 
boostinvariant hydrodynamics according to their eigentime). While such a local freeze- 
out criterion would be a more realistic scenario, it is very difficult to realize, especially in a 
fully three-dimensional calculation. 

Due to these difficulties and in order not to introduce additional free parameters (like the 
freeze-out temperature and / or density), we prefer to discuss only fluid instead of particle 
distributions. Our results (e.g. the dN/dY-spectra, cf. section |4.2| ) are therefore compared 
to the one-fluid limit, not directly to experimental data. Also due to this reason, we show 
the time evolution of the observables under consideration at fixed center-of-mass time (not 
at fixed proper time). 
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4.1 Density evolution 



In section |2.1| we claimed that in a one-fluid model much higher densities are reached than 
in multi-fluid models. Figure [l] shows the average baryon density of the projectile nucleus 
as a function of time. The averaging is defined as 

Jd: 3 xRp(x)np(x) 
{Up) ~ JVxR P (x) ■ (29) 



where Rp denotes the projectile baryon density in the computational frame, cf. section 2.2 



Thus it is a measure for the compression of the projectile only. In the one-fluid limit the 
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Figure 1: Average baryon density of the projectile. 



maximum projectile density is almost a factor of two higher and is reached earlier than in 
the three-fluid model. 

The dotted line shows a calculation without fireball production (erg = 0), like in the one- fluid 
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limit, but the unification was performed after 2.5/m/c. The evolution first follows more or 
less the three-fluid curve but yields slightly higher densities. It was already pointed out in 
f30fl that neglecting the energy loss due to particle production results in higher excitations in 
multi-fluid models. At 2.5/m/c, when instantaneous unification like in the one-fluid limit is 
applied, the curve jumps towards the one-fluid limit but does not reach the full height since 
the nuclei already penetrated partially. 

The reason for the lower compression of the baryons in the three-fluid model is the finite 
stopping length of nucleons in nuclear matter. This is illustrated in Fig. |2| which shows the 
energy density profile of the projectile fluid along the beam axis. In the one-fluid limit the 
shock wave which moves outwards can be clearly seen, leaving behind compressed matter 
with energy density e ~ 20eo- Also, no flow through the plane z = exists (corresponding to 
instantaneous, complete stopping of the projectile). We checked that the shock- wave velocity 
and the energy density of the shocked matter in this three-dimensional calculation agrees 
satisfactorily well with the one-dimensional shock model. In contrast, in the three-fluid 
model there is no discontinuity in the baryon or energy density profile which separates the 
central region from the incoming ground-state nuclear matter. Due to their non-vanishing 
mean-free path some projectile nucleons do pass the (z = 0)-plane and the shock front is 
broadened considerably. As a consequence, the maximum energy density of the projectile 
matter does not exceed 10eo- 

4.2 Baryon stopping 

Figures || and [| compare the time evolution of the rapidity distributions of the nucleonic 
fluids 1,2 in the one-fluid limit and the three- fluid model, respectively. We again stress that 
both figures show fluid distributions, i.e. the particle momenta are not thermally smeared. 
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Figure 2: Energy density profile (in units of e ~ 0.15 GeV/fm 3 ) of the projectile fluid along 
the beam axis (in the plane x — y — 0). 
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The time steps in the two figures are not equal, since the compression reaches its maximum 
faster in the one- fluid limit, cf. Fig. |l|. 



The major difference between one- and multi-fluid models as claimed in section 



is 



becoming obvious: At 1.2 fm/c the fluid rapidity distribution in the one-fluid limit (Fig. [| 
upper left panel) shows clearly that participant matter stops immediately at midrapidity 
when the nuclei touch, with nothing in between midrapidity and the projectile or target 
region, respectively. Comparing the rapidity distributions at 4.8 fm/c (Fig. middle right 
panel, and Fig. f|, upper right panel), in the one-fluid limit the colliding baryons are piling 
up at midrapidity while in the three-fluid model they are shifting smoothly towards midra- 
pidity. Nevertheless, looking at the later stages, i.e. 6.0 fm/c, even in the three-fluid model 
the colliding nuclei become more or less stopped within about one unit in rapidity before 
reaccelarating. Thus, the coupling sourceterms, parametrized by free binary nucleon-nucleon 
collisions, are strong enough to yield a sufficient stopping. 

A benefit of our unification procedure is that we can track the individual nucleonic fluids 
even after unification since they are formally propagated on separate grids. We observe that 
projectile and target are changing position in the dN/dY in the one-fluid limit, while they 
keep the sign in rapidity in the three-fluid model, i.e., that projectile and target are reflected 
in the one-fluid limit, while heading on in their initial direction in the three-fluid model. 

In the late stage of the reaction expansion is becoming the most important process which 
is reaccelerating the fluids again in both the one-fluid limit and the three-fluid model. 



4.3 Equilibration 

For an equilibrated system the pressure is isotropic. In contrast, while equilibration is not yet 
established, the transverse and longitudinal pressures (as defined in kinetic theory ||35|| ) will 
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differ from each other. In our approach, this reflects in a deviation of the sum of the partial 
pressures (of projectile and target fluids) from the equilibrium pressure. This is investigated 
in Fig. ^ by calculating the ratio p± + p 2 to p equ u that a unification of the nonequilibrated 
fluids would yieldQ. The same analysis is done for the rest frame baryon densities, 

^3^, pi(x)+p 2 (x) r^3 x ni (x)+n 2 {x) 

f Pl + P2 ^ _ J Pequil(x) + _ J n equ g(x) {^S) 

Pequil fd^X ^equil fcfix 

The volume integration runs in one case only over a thin central layer perpendicular to the 
beam axis, so that most of the matter traverses this plane during the reaction and in the 
second case over the whole overlap region of fluids 1 and 2. The ratio is by definition equal to 
unity in the case of unified fluids, for in that case p\ and p 2 just denote the partial pressures. 
One observes that kinetic equilibrium is established only after a time t^ M ~ (5 — 6) fm/c ~ 



2Rau/icm- A similar equilibration time scale was also extracted from a cascade model [36fl. 



Thus, observables which are sensitive to this early stage of the reaction |37| (as e.g. the 
nucleon flow discussed in the next section) will deviate considerably from their values in the 
one-fluid limit. One also sees (cf. dashed curve) that the merging of projectile and target 
nucleons starts at tcu ~ 5 fm/c, when the collision terms F^ 2 have almost established 
kinetic equilibrium. 

It seems that the matter in the central plane equilibrates faster than the whole overlap 
region. However one has to recall that matter is passing through this thin layer, so that the 
ratios evaluated at different times refer to different parts of matter. Moreover, the respective 
matter corresponds to matter trailing the fronts of the nuclei. On the other hand, for the 
overlap region more or less the same matter remains within the volume of integration at all 
times. 



1 Pequil and rieguii are calculated from eqs. (0) by inserting E = E\ + E 2 , M = M\ + M 2 and R= R\+R 2 
for the left-hand side. 
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Figure 5: Evolution of sum of nonequilibrium pressures and rest frame baryon densities 
compared to the corresponding equilibrium values (ratios). 
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4.4 Nucleon flow 

As a measure for the in-plane transverse flow we investigate the mean transverse momentum 
per nucleon (p x /N) as a function of rapidity y. We define (p x /N) as follows: 

, , , x f„d 3 a; R(x) mjv u x (x) , 

= - sXro) (31) 

Here u x = jv x denotes the component of the local 4-velocity field in the direction of the 
impact parameter and ttin the nucleon rest mass. We assume m^u x to be a measure for 
the component of the nucleonic momentum in x-direction without thermal smearing. The 
volume integration is performed over all fluid elements (projectile and target) around a given 
rapidity y. Figures |] and |7] show the time evolution of the (p^/A^-distributions for one-fluid 
limit and three-fluid model, respectively. 

During the reaction pressure builds up and "spectator" matter is deflected. Therefore, 
(Px/N) is increasing. The one-fluid limit shows higher maximum (p x /N) than the three-fluid 



model, i.e. produces more sideward flow, as claimed in section |2.1| . In the later stages the 
flow is still ~ 280MeV/c in the one-fluid limit, which is higher than the final value of the 
three-fluid model of ~ 120MeV/c. The flow for the three-fluid model does not change any 
more from about 6.0 fm/c on and thus is a reasonable measure for the final amount of flow 
after freeze-out. 

In [10] an excitation function of the directed flow (p x hr ) has been investigated for a 



hadron gas EoS in comparison to an EoS with a phase transition to QGP. This one-fluid 

of ) 



model predicted a minimum of the directed flow \p d x ) around the phase transition region. 



Px) was defined as [10 



p * r ) = l I dy ^w\ iiPx/NHy)) • (32) 



This is basically the weighted mean of the function {p x /N) (y), eq. (pl|). 
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Figure 6: The evolution of (p x /N) in the one-fluid limit. 
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Figure 7: The evolution of (p x /N) in the three-fluid model. 
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For our EoS QT7p, (pf r ) = 120 MeV/c in the one-fluid limit, whereas (pf r ) = 40 MeV/c 
in the three-fluid model. Even with our crude EoS without a phase transition ( |i~7|) the di- 
rected nucleon flow is as low as in the one-fluid model with a phase transition to QGP [ |T(| . 
Obviously, the effect of the nonequilibrium effects is of the same order of magnitude as 
the softening of the EoS due to a (first order) phase transition. Some implications of ki- 
netic nonequilibrium in the early stage of ultrarelativistic heavy-ion collisions were already 
sketched in [07 . 



5 Summary, conclusion and outlook 

In this paper we presented a three-fluid hydrodynamical model which allows to account for 
nonequilibrium effects between target, projectile, and produced particles during the early 
stage of the collision. We discussed that due to the nonvanishing thermalization time scale, 
this model yields a 

1. lower energy density and temperature, 

2. lower pressure, 

3. less baryonic compression, 

4. and a different velocity distribution 

of the nucleons at early times as compared to the one-fluid hydrodynamical model (which as- 
sumes instantaneous local thermalization between projectile, target, and produced particles). 
As a consequence, we expect that the 

1. directed nucleon flow, 
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2. entropy production (pion multiplicity), 



3. thermal electromagnetic radiation |[21|| , 



4. lifetime of the hot and dense central region 

differ considerably in the three-fluid as compared to the one-fluid model. These results 
suggest that the predictions of the one-fluid model (as summarized in the introduction) have 
to be modified by taking nonequilibrium effects into account. 

We especially point out that in three-fluid model the projectile and target nucleons reach 
kinetic equilibrium only after t e J M rs 2R/jcm- Thus the directed nucleon flow (emerging 
at early times) is smaller and in particular much less sensitive to the EoS as compared to 
the one-fluid limit. This is due to the fact that the flow of the nucleons of fluid i is caused 
mainly by the partial pressure pi and not by the full p eq (see Fig. [|). 

In the future more realistic equations of state for the projectile and target fluids (in- 
corporating a phase transition to a quark-gluon plasma at high energy densities) have to 
be implemented and the excitation functions for the pion multiplicity, the mean transverse 
momentum, and the directed nucleon flow have to be reinvestigated. 
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